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ABSTRACT 

In recent years, analyses of eclipsing binary systems have unveiled differences between the observed 
fundamental properties of low-mass stars and those predicted by stellar structure models. Particularly, 
radius and effective temperatures computed from models are ~ 5 — 10% lower and ^ 3 — 5% higher 
than observed, respectively. These discrepancies have been attributed to different factors, notably 
to the high levels of magnetic activity present on these stars. In this paper, we test the effect of 
magnetic activity both on models and on the observational analysis of eclipsing binaries using a 
sample of such systems with accurate fundamental properties. Regarding stellar models, we have 
found that unrealistically high spot coverages need to be assumed to reproduce the observations. 
Tests considering metallicity effects and missing opacities on models indicate that these are not able 
to explain the radius discrepancies observed. With respect to the observations, we have tested the 
effect of several spot distributions on the light curve analysis. Our results show that spots cause 
systematic deviations on the stellar radii derived from light curve analysis when distributed mainly 
over the stellar poles. Assuming the existence of polar spots, overall agreement between models and 
observations is reached when ~ 35% spot coverage is considered on stellar models. Such spot coverage 
induces a systematic deviation in the radius determination from the light curve analysis of ~ 3% and 
is also compatible with the modulations observed on the light curves of these systems. Finally, we 
have found that the effect of activity or rotation on convective transport in partially radiative stars 
may also contribute to explain the differences seen in some of the systems with shorter orbital periods. 

Subject headings: binaries: eclipsing - binaries: spectroscopic ~ stars: late-type - stars: activity - 
stars: spots - stars: evolution 



1. INTRODUCTION 

I> I Analyses of double-lined eclipsing binary stars (here- 
^D , after EBs) have provided fundamental properties of stars 
CN on the low-mass end of the Main Sequence. Such EB 
r ' systems have revealed to be the best candidates to test 
^'j , stellar structure models since they can potentially yield 
l/~j ■ masses an d radii for M-type sta rs with uncertainties be- 
^^ ' low ~ 1% ([Morales et al.l[2009al ). However, results com- 
<^ ' ing from a variety of independent studies indicate that 
T— I . differences exist between the observed properties and 

ILJ ■ those predicted by models; for the same mass (M), mod- 
. !^ ] els predict ~ 5 — 10% smaller radius (i?) and ~ 3 — 5% 
S^ i larger effective temperatures (Tog) than observed for 

;_( ' stars below 1 M0. In contrast, the stellar luminosity 

Ci I for these binary components is in agreement with those 
of sin gle stars and also with theoretical models (see Ribas 
I2OO6I for a review). 

Several authors have argued that activity could 
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be the cause of discrepancies found between mod- 
els a. n d observations of these syste i ns ([Torres &: RibasI 
200I lL6pe^-Morales k Ribasl [20051: [Torres et all 120061: 



Lopez-MoralesI I2007t iRibas et all l2008| ) . Further sup- 
port to this hypothesis is lent by similar radius dis- 
crepancies found when comparing single magnetically 
active late-type sta rs with their inactive counterparts 
([Morales et al.|[2008[ ). and by the results o f interferomet- 
ric o bservations of single inactive stars ([Demorv et al.l 
|2009| ). Until now, all the known low- mass EBs are close 
binaries with periods well belo w 10 days, th us they are 
expected to be synchronized ([Mazehl [2008[ ). i.e., with 
their rotation tidally locked to the orbital motion, and 
therefore they are fast rotators (urot ^ 10 km s~^). Such 
fast rotation in presence of magnetic fields could trig- 
ger high levels of activity that may modify the struc- 
ture of these stars in some respects. Observationally, the 
high activity of these stars is usually confirmed by sat- 
urated X-ray emission {\ogLx/Lboi ~ —3), strong chro- 
mospheric emission in the Ha Balmer line, fiares, and 
modulations on the light curves caused by the presence 
of starspots on the stellar surface. 

Over the past decades, substantial progress towards 
the understanding of the structure and evolution of low- 
mass stars has been made, also including activity ef- 
fects. iChabrier et al.l ([2007D showed that the introduc- 
tion of rotation and /or ni agnetic field effects on the mod- 
els of lBaraffe et al.l |l998[) could explain the discrepancies 
between the observed and theoretically predicted mass- 
radius relationship of EBs. Specifically, they presented 
two scenarios considering: (1) that the effect of mag- 
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netic field and rotation alter the efficiency of convective 
energy transport, which can be modelled by setting the 
mixing length parameter (a from Mixing Length Theory) 
to lower values than those used for solar models; and (2) 
that the stellar magnetic activity present on these objects 
can be associated with the appearance of dark spots cov- 
ering the radiative surface, which can be modelled by 
assuming a new stellar luminosity, L ex (1 — /3) R^T^ff, 

where R and Tcff are the modified stellar radius and ef- 
fective temperature, respectively, and /3 is the factor of 
spots blocking the outgoing luminosity due to their lower 
temperature (see Section [3] for further details on this pa- 
rameter). 

The results of iChabrier et alj ()2007D show that both 
these scenarios predict larger radius than standard stel- 
lar models, but while the effect of spots is significant over 
the entire low- mass domain, the effect of convection is rel- 
atively smaller for fully convective stars (M < OAMq) 
because convection is nearly adiabatic and changing the 
mixing length parameter has a modest impact on their 
stellar structure. Since radius discrepancies have been 
reported over a wide range of masses, this may be indi- 
cating that either the contribution of spots or a combina- 
tion of spots and activity effect on convective transport 
are needed to reconcile models with observations. 

New accurate fundamental prop erties of low-mass EBs , 
and especially the case of CM Dra (jMorales et al.|[20Q9aD . 
providing precise characterization of two fully convec- 
tive stars, makes it worthwhile to carry out a new thor- 
ough analysis of the agreement between observati on and 
theory using the scenario suggested in Chabric r et al.l 
(|2007t) . In Section [5] we briefly describe the discrepan- 
cies between stellar structure models and current obser- 
vations of low-mass stars and we select the sample of 
best known EB systems to constrain stellar models. In 
Section [3] we analyze the effect of the photospheric spots 
on EB light curves and on the parameters fitted to these 
curves. Critical comparison with models including ac- 
tivity effects is discussed in Section |4l and finally, the 
conclusions are presented in Section [5] 

2. LOW-MASS EBS COMPARED TO MODELS 

The properties describing the evolutionary phase of a 
star are its mass, chemical composition and age. Given 
these three properties, radius and effective temperature 
are determined from stellar evolutionary models. There- 
fore, models could be tested comparing the theoretical 
M — R and M—Tes relationships with observational val- 
ues. In Figure [J such comp arisons are shown for the 
models of lBaraffe et al.l ()1998[ ) and the full sample of low- 
mass stars with known physical properties. This figure 
includes properties derived from double-lined and single- 
lined EBs, interferometric radius measurements and also 
from analysis of M-type stars in post-common envelope 
binaries, although some of these measurements are not 
independent of theoretical models or calibrations. In 
post-common envelope components, the systems may not 
be reliable representations of the stars described by mod- 
els due to interactions suffered during the common enve- 
lope phase. 

In Figure [TJ a clear difference between stars below and 
above ~ 0.35 M0 is apparent. On the lower end, radii 
are much closer to theoretical models and less scattered. 



while on the upper end, the scatter is larger and there 
seems to be a larger difference with models as well. This 
was already noticed in previous analyses (jRibasI [2006D 
and it may be consistent with the boundary between 
fully convective stars {M < 0.35 M0) and partially ra- 
diative stars. Considering these two regimes, differences 
between all known measured radii of low-mass stars and 
the 1 Gyr solar metallicity model isochrone are 5.7% and 
8.1% for fully convective and partially radiative stars, 
respectively. The dispersion of these differences is 3.1% 
and 12.2%, respectively, thus illustrating the scatter dif- 
ference mentioned above. In the case of effective temper- 
atures the differences are —1.9% and —1.6%, respectively, 
for this sample. 

To test stellar structure models, masses and radii of 
stars need to b e known with a c curacies down to th e few 
percent level (Andersen 19911; iTorres et al1l2009f ). As 
mentioned earlier, EB systems provide the best opportu- 
nity to measure these fundamental stellar properties with 
such accuracy, as clearly visible in Figure [1] Thus, for 
our analysis we have selected all known low-mass stars 
with accuracies better than 3%. As reported recently 
bylTforres et al. (2009), special attention should be paid 
both on the precision and the accuracy of the measures, 
and so we have excluded from our sample systems with 
error estimations not clearly described. Furthermore, we 
have only selected main sequence systems with ages be- 
low 5 Gyr and masses below 0.8 M0 to restrict radius dif- 
ferences due to evolution below 2% and also to minimize 
the effect of mixing length differences between low-mass 
and solar mass models. The mechanical and thermal 
properties for each component of the selected systems 
are listed in Table [1] as well as the age and metallicity 
when known. Comparison of this sample with models is 
shown on the right panels of Figure [TJ For the subsample 
of Table [H the comparison with the 1 Gyr solar metallic- 
ity model gives discrepancies with observations of 5.8% 
and —6.1% in radii and effective temperatures, respec- 
tively, for fully convective stars (CM Dra), and 7.8% and 
—3.0%, respectively, for partially radiative stars. 

Besides magnetic activity, severa l model parameters 
such as metallicity or opacities (|Berger et al.l 120061 : 
iCasagrande et al.l |2008() have been proposed to explain 
the discrepancies between theory and observations. Our 
calculations indicate that metallicity produces differences 
between ~ 4% and ~ 20% on the radius and effective 
temperature values predicted by models, respectively, 
when comparing [m/H] = 0.0 and [m/H] ~ —1.0 models 
(solid and dashed line, respectively, in Figure |T|). How- 
ever, none of the systems is known to have such low 
metallicity, except for the case of CM Dra, for which its 
metallicity is not well estabhshed, and NGC2204-S892, 
for which there is no determination. Besides, radius dif- 
ferences are lower than needed to reproduce observations 
even with metallicities as low as [m/H] = —1.0. There- 
fore, metallicity effects are expected to be not significant 
in our sample. In the case of effective temperature com- 
parisons, where the effect of metallicity is more impor- 
tant, EB light curves only provide accurate values of the 
temperature ratio for the components. The individual 
effective temperatures are usually computed using dif- 
ferent spectrophotometric calibrations for each system, 
thus introducing the potential for systematic errors. Fur- 
thermore, some EB systems may have peculiarities, such 
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Fig. 1. — Different sets of theoretical M — R and M — T^ff relationships compared with observational values of low-mass stars. Left: 
Low-mass stars in double-lined EB systems (filled circles), single-lined EB systems (open circles), M stars in post-common envelope binaries 
(grey squares) and late-type stars with radius measured interferometrically (open triangles). Right: Same plots showing the E Bs included 
in our analysis (see text and Table [1} . Theoretical models with different sets of parameters as labelled are from IBaraffe et al.| |) 1998| ) . 



as CU Cnc, which is supposed to have a circumbinary 
dust disk that may affect the temperature determination. 
Tlius, to investigate the agreement between models and 
observations we will primarily make use of the M — R 
relationship, which are both fundamentally determined 
quantities. 

Recently, the possibility of underestimated opacities 
has been raised in the context of the Standard Solar 
Model. An increase of opacities near the bottom of the 
convective zone has been shown to reconcile helioseis- 
mic data of t he Sun with models u sing t he new chemical 
composition (|Bahcall et all 12004 120051 ). These results 
and others found in the literature invoking problems with 
the atmospheric opacities as a source of the discrepancies 
present in observed radii of EB compared with theoret- 
ical models prompted us to perform different tests to 
analyze the effect of atmospheric opacity on the result- 
ing stellar radius. Using the standard Lyon group stel- 
lar evolutionary models (jBaraffe et al.lfl998l ). an overall 
massive increase by a factor of '^ 10 in the opacity affects 
the radius by just 4% (and even this difference is smaller 
be low the fully convec t ive bo undary), as already shown 
bv lChabrier fc Baraffd (|T997l ). To obtain a ~ 10% larger 
radius the opacities of eclipsing binaries would have to 
be increased to unrealistically high values. Missing el- 



emental opacities do not seem to be enough to explain 
the '^ 5 — 10% radius disagreement between models and 
observations. Moreover, both metallicity and opacity ef- 
fects should also be present on single stars, which are 
known to be well described by models, as recently con- 
firmed by iDemorv et al.l (|2009l ) from the comparison of 
interferometric radius measurements of inactive low- mass 
stars with model predictions. Further, the authors did 
not find a correlation between metallicity and radius dif- 
ferences thus not confirming the results from Berger et 
al. (2006). Further, similar radius differences as those 
reported here are known between single m agnetically ac- 
tive s tars and their inactive counterparts ([Morales et al.l 
l2008f l. thus strengthening the activity hypothesis. 

3. DETERMINATION OF PHYSICAL PROPERTIES OF EBS 

The conclusions of the previous section reinforce mag- 
netic activity as the most reasonable explanation for the 
discrepancies found between stellar models and obser- 
vations of low-mass stars. Before embarking on a care- 
ful fine-tuning of model parameters to reach agreement 
with observations it is worth evaluating whether our ob- 
served eclipsing binary parameters could be affected by 
the strong magnetic activity of the components. Time- 
variable brightness changes are a common signature of 
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TABLE 1 

Properties of observed EBs included in our sample. 





P 


M 


R 


Teff 




Age 


f rot, sync 




EB 


(days) 


(Mo) 


(R©) 


(K) 


[M/H] 


Gyr 


(kms-i) 


Ref.^' 


V818 Tau B 


5.61 


0.7605±0.0062 


0.768±0.010 


4220±150 


0.13 


0.6 


6.95±0.09 


1 


IM Vir B 


1.31 


0.6644±0.0048 


0.681±0.013 


4250±130 


-0.28 


2.4 


26.12±0.50 


2 


NGC2204-S892 A 


0.45 


0.733±0.005 


0.720±0.010 


4200±100 


- 


- 


80.9±1.1 


3 


NGC2204-S892 B 




0.662±0.005 


0.680±0.020 


3940±110 






76.4±2.2 




GU Boo A 


0.49 


0.6101±0.0064 


0.627±0.016 


3920±130 


-0 


1.0 


64.7±1.6 


4 


GU Boo B 




0.5995±0.0064 


0.624±0.020 


3810±130 






64.4±1.6 




YY Gem A & B . 


0.81 


0.5992±0.0047 


0.6194±0.0057 


3820±100 


-0 


0.4 


38.88±0.36 


4 


GU Cnc A'' 


2.77 


0.4349±0.0012 


0.4323±0.0055 


3160±150 


-0 


0.32 


8.02±0.10 


4 


GU Gnc B'' 




0.3992±0.0009 


0.3916±0.0094 


3125±150 






7.30±0.17 




GM Dra A 


1.27 


0.2310±0.0009 


0.2534±0.0019 


3130±70 


-1/-0.6 


4.1 


10.02±0.08 


4 


GM Dra B 




0.2141±0.0008 


0.2398±0.0018 


3120±70 






9.51±0.07 





'^ References: 1. [Torres fc Rlbai 1(2002 ): 2. [Morales et"all H2009bl 1: 3. IRozvczka eTaP l[200gf ): 4. [Torres et all l[2009f ) 
^ T(,ff could be underestimated due to the presence of circumbinary dust. 
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the presence of activity on low-mas stars both on single 
and eclipsing binaries. In the latter case, the modula- 
tion present on the light curves is combined with the 
eclipsing variability and, therefore, starspots must be 
taken into account to derive accurate stellar fundamen- 
tal properties. In this section we study the possibility 
that starspots could be responsible for systematic effects 
inherent to the light curve and radial velocity analyses 
in the case of the very active low-mass eclipsing bina- 
ries. We analyze whether the fundamental properties 
determined from the classical modelling could be biased 
because of the presence of starspots in different geome- 
tries. 

Absolute physical properties of EBs are obtained 
from the combination of the results of light and ra- 
dial velocity curve analyses with modelling codes such 
as the widely-used Wil son-Devinney (hereafter WD; 
IWilson fc Devinnevlll971[ ). Light curves yield properties 
such as the inclination of the orbit (z), the ratio of effec- 
tive temperatures {Tcff2/Tosi), the luminosity ratio in 
the light curve bandpass (L2/-L1) and the surface poten- 
tials (i7i and J72), which are related to the radii relative 
to the semimajor axis (ri and r2). Radial velocity curves 
provide values for the mass ratio (g) and projected semi- 
major axis (asini). The combination of the two mod- 
elling procedures yields the absolute properties of the 
components independently of models or distance calibra- 
tions. Radiative parameters, such as limb darkening and 
gravity darkening, relevant for light curve analyses, are 
usually taken from theory. 

The presence of spots on the surface of the components 
of an EB system cause perturbations on both the light 
and radial velocity curves. The most prominent effect 
is on light curves, in which modulations in the out-of- 
eclipse phases appear due to the transit and occultation 
of spots according to the orbital and spin motions. To ap- 
propriately derive physical properties of the components, 
spot effect must be taken into account when analysing 
these curves. The WD code introduces spots on the mod- 
elling assuming that they are circular and that they have 
an uniform temperature ratio with respect to the photo- 
sphere. With this simple model, properties of spots, such 
as their location in co-latitude {9s , measured from North 
pole) and longitude {(jjs), the angular size (rg) and the 

temperature contrast with the photosphere, i.e., j; '" , 



where Tcs and Tcs^s are the effective temperatures of the 
photosphere and the spots, respectively, can be fitted. 
The values of 0s and Tg are well constrained by the central 
phase of the modulation and its duration and amplitude, 
respectively, however the amplitude is also dependent on 

0s and J ;° , so these parameters are strongly correlated 
with the size of spots. This model can reproduce the 
presence of complex and irregular spot groups by equiv- 
alent circular spots with an average temperature. 

The total luminosity of a spotted star is the addition of 
the contribution of the spots at an effective temperature 
Toff'^s and the immaculate surface at Tog. If we consider 
S and Ss to be the surface of the star and that covered 
by spots, respectively, the total luminosity is given by: 



L^iS^ Ss) <jZ 



cff 



SsCrTos , 



(1) 



Co mparing this equation with the formalism introduced 
bv lChabrier et al.l (|2007[ ). it can be concluded that the (3 
parameter is given by: 



^ S 



1- 



'-eS. 



T, 



off 



(2) 



This equation indicates that the /3 parameter is related 
with the properties of spots. For completely dark spots 
(T'cff,s = 0), the term in brackets in Equation[2]is exactly 
one, therefore, /3 is a measure of the fraction of stell ar sur- 
face c overed by dark spots as defined by Chabric r et al.l 
(|2007D . However, in the realistic case that spots are not 
completely dark the term in brackets is less than one 
and /3 is lower than the fraction of spotted surface. For 
a given T^s and an estimation of the contrast of spots. 
Equation [5] provides the ratio of surface covered by spots 
when P is known, or vice versa. 

Spot modelling results have b een published for sev- 
eral of the best-known EBs (iTorres fc Ribai 120021: 



Riba! 



2009 



gOOaiLopez-Morales fc Ribasll2005l: [Morales et all 



3) reporting both cold an hot spots with radii be- 
tween ^ 9 deg and ^ 90 deg and with temperature fac- 
tors (Toff s/Tcff) down to 0.86. Hot spots can be in- 
terpreted as photosphere regions surrounded by large 
cool spots. Considering the spot parameters given in 
the above references for the cases of IM Vir, GU Boo, 
YY Gem, CU Cue, and CM Dra, we find P values of 
up to 0.1 according to Equation [2] This value is much 
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smaller than th e rang e between 0.3 and 0.5 suggested by 
iChabrier et al.l ()2007l ). It must be mentioned, though, 
that since the relevant measures in EB light curve anal- 
yses are differential magnitudes, the photometric varia- 
tions used to derive the spot parameters are not sensitive 
to the total surface covered by spots but to the contrast 
between areas with different effective temperatures. For 
instance, an evenly spotted star would not show signifi- 
cant light curve variations. It is thus possible that EBs 
could be more heavily spotted than simple photometric 
variations indicate and permit /3 values higher than 0.1. 
A check of the photometric variations that could produce 
such higher values of /? and their effect on derived param- 
eters from fits are described in the following sections. 

3.1. Modulation on light curves 

For M-type stars, an assumption of roff"^3500 K and 
spots about ^ 500 K cooler than the photosphere and /? 
values between 0.3-0.5 corresponds to a fractional sur- 
face area of 65-100% covered by spots (as shown by 
Equation [2]) . To test the consistency between the the- 
oretical modulations expected from different /? values 
and the spot modulations observed on EB light curves, 
we simulated EB systems with components randomly 
spotted and used the WD code to generate their light 
curves. Such theoretical light curves were subsequently 
compared with those observed in EBs. 

We have developed a code to randomly place spots 
on the surface of stars. We assumed a uniform longi- 
tude dist ribution and t r ied d ifferent distributions over 
latitude. iGranzer et al.l ()2000[) calculate the probability 
of spot appearance for different latitudes as a function 
of spin period and mass. Their results show that for 
a star like the components of GU Boo (see Table [1]), 
spots are formed in a band from 25 deg to 55 deg of 
co-latitude. Both bands (one per hemisphere) repre- 
sent about one third of the stellar surface, which is 
not enough to simulate spot coverages with /? greater 
than 0.15 (assuming spots 500 K cooler than the pho- 
tosphere). Doppler tomography analyses have revealed 
the prominent exis t ence of polar spots on a ctive stars 
(iJeffers et al.l 120071: I^Vashuettl fc StrassmeieU 2001 : see 
also iStrassmeieii 120091 and references therein) . Thus , 
we decided to extend the distribution to all stellar lat- 
itudes. To use simi lar distributions as those given in 
IGranzer et al.l (|2000l ) we considered a linear probability 
from the pole to 40 deg and from 45 deg to the equa- 
tor with a plateau between 45 deg and 40 deg (hereafter 
Distribution 1). This roughly mimics distributions of 
0.6 Mq rapidly-rotating stars. For comparison we consid- 
ered also a bilinear distribution from the pole to 70 deg 
with a peak a t 25 deg. Thi s is similar to the 0.4 M© 
distribution in IGranzer et al.l (|2000) for rapidly rotating 
stars but extended to cover /? factors up to 0.3 (hereafter 
Distribution 2). Finally, we considered a completely uni- 
form distribution. 

These distributions are shown in Figure [2] and it illus- 
trates that Distribution 2 is more concentrated towards 
the poles. Additionally, to have reasonable computing 
time using the WD code, we limited the number of spots 
by imposing that the centers of two spots could not be 
closer than half their radius, i.e., each spot contributes 
at least ^ 30% of its surface to the total spot coverage. 

The simulation procedure was executed as follows. We 



3.2 
2.8 
2.4 

d 
"^ 1.6 

p 
3- 
1.2 

0.8 

0.4 



n l-i- -'i 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

cos(O) 

Fig. 2. — Probability density functions over cos 9 used to simu- 
late spotted stars. The equator corresponds to 6 = 90 deg. The 
uniform distribution (solid line), Distribution 1 (dot-dashed line) 
and Distribution 2 (dashed line) are described in the text. 



considered spots with a temperature contrast with the 
photosphere of 0.85, even on overlapping areas (i.e., the 
temperature contrast in the area where two spots over- 
lap is also 0.85). This value constrains /3 between and 
~ 0.5, the latter meaning a completely spotted stellar 
surface. In the case of Distribution 2 the upper value of 
(3 is ^ 0.35. We chose 10 deg as spot radius in line with 
large spot groups measured on the Sun. Then, for each 
latitude distribution we simulated 25 light curves with 
spots randomly distributed over the surface of the com- 
ponents for different values of f3. From these simulated 
light curves we measured the peak-to-peak variations on 
the out-of-eclipse phases and computed the mean value 
of the 25 light curves. This procedure was done using 
the physical properties of GU Boo as template and using 
R band light curves. However, since the effect of abso- 
lute masses or mass ratio on the light curves is negligible, 
this should be representative of the entire low-mass star 
domain. 

The values of the peak-to-peak modulations found with 
these simulations were later compared with the real val- 
ues for the best-known EBs. Figure [3] shows the peak- 
to-peak modulations resulting from simulations for each 
distribution compared with a few observational values 
as reference. Large modulations are better explained 
by mild spot coverages assuming a uniform distribution 
or Distribution 1, while low spot signals are better re- 
produced by Distribution 2. However, the amplitude 
of these modulations depends on the size and thermal 
properties of the spots as shown on Figure |4l It is clear 
that the observed modulations could be explained with 
any of the distributions when considering the different 
spot sizes and the intrinsic dispersion of the modulation. 
Therefore, low amplitudes such as those of YY Gem and 
CM Dra could be recovered with a mild spot coverage 
with small size or low temperature contrast spots. Be- 
sides, the variability of the spot properties should also be 
taken into account since it could induce different modu- 
lations on different epochs as has been reported for the 
case of V band light curves of YY Gem, showing am- 
litudes of 0.09 mag a nd 0.055 mag in different seasons 
Torres fc Ribai[2002l ). 

Several other distributions, such as linear, quadratic 
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Fig. 3. — Amplitude of the modulations produced by different 
spot scenarios using different spot distributions on R band light 
curves. Each point is the mean value of 25 realizations and the 
shaded area shows the standard deviation of these modulations. 
The values reported in the literature for CM Dra (dashed line), 
YY Gem (dot-dashed line) and GU Boo (long-dashed line) are 
shown for comparison. 

and square root from equator to poles were tested yield- 
ing results between the uniform and Distribution 1 cases. 

3.2. Systematics on light curves due to spots 

The simulations described above, were also used to 
investigate systematic effects of spots on fits to light 
curves. A thorough statistical analysis of the system- 
atic deviations induced by the presence of spots on these 
parameters would require fitting all the simulated light 
curves by treating them as real observations. Typical 
light curve analyses involve testing fits for different pa- 
rameters, especially when spots are present, since differ- 
ent spot configurations should be tested. Simultaneous 
fits of light curve in several bands are also preferable 




Fig. 4. — Same as Figure[3]for the uniform distribution comparing 
the cases of simulations using spots with radius of 15 deg (light 
shade) and 5 deg (dark shade). 

in order to avoid correlations between parameters such 
as limb darkening and radius. Besides, spot parameters 
are degenerate and thus the inversion of a light curve 
to recover the input parameters is an ill-posed problem. 
This implies carrying out a slow procedure consisting in 
adding spots and evaluating the relevance of the mod- 
ulation to reach the best compromise between stability 
and the quality of the fit to the observations. For these 
reasons, fitting all the simulated curves is prohibitively 
time consuming. We performed fits to three simulated 
light curves for each distribution and for each of the /3 
values 0.1, 0.2, 0.3 and 0.4 (33 light curves in total) with 
the aim of inferring a general trend. Besides, we focused 
our attention on the i?-band to reduce the computing 
time. This band is commonly used in the photometric 
follow-up of low-mass eclipsing binaries and it is thus 
representative of light curve analyses. A possible bias 
resulting from the correlation between radius and limb 
darkening coefficients should be negligible because the 
coef ficients are co mputed from the same theoretical ta- 
bles (jClaret ' 2000) both for simulations and fits. For each 
of the f3 values mentioned, we chose as representative 
spotted light curves three realizations with modulation 
amplitudes representative of the average. 

We added random Gaussian noise of 1% of the flux 
to the simulations so that the typical scatter of light 
curves is reproduced. We subsequently performed fits 
with the WD code. Light curve relevant parameters such 
as the inclination, temperature ratio, light ratio and sur- 
face potentials were fitted along with spot properties. 
Since some of the spot parameters are highly correlated 
we started assuming spots at 45 deg latitude and 10% 
cooler or hotter than the photosphere. The free param- 
eters were the longitude {(j)s) and size (rg) of spots and, 
in a second step, the latitude {9s) and temperature con- 
trast (Teff g/Teff). As the usual practice, we included in 
the model as many spots as needed to obtain a realis- 
tic fit and we tested the location of spots on primary 
component, secondary or both at the same time. All 
cases yielded satisfactory fits with 1-3 spots on different 
components. We found best fits with Tes^g/Tes ratios of 
~ 0.9, thus indicating that spot groups could be mod- 
eled with a smoother temperature contrast than real. A 
phase shift was also set as a free parameter to account 
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Fig. 5. — Two representations (top panel for Distribution 1 and bottom panel for Distribution 2) of simulated spotted EB systems 
(primary component on the left and secondary on the right) with the resulting light curves. Star surfaces are represented in a MoUweide 
projection and the center of each one corresponds to deg of longitude and 90 dcg of co-latitude (equator). Best fits to the simulated light 
curves together with the O-C residuals are shown. Both cases were simulated with /3 = 0.2. The reference level of the unspotted light is 
shown for comparison (dot-dashed line). 



for phase corrections due to spots. 

Figure [5] depicts two examples of the fits together with 
the configuration of the simulated spots on the stars, 
and Figure [6] shows all the 33 fits performed. In Figure [7] 
we illustrate the effect of spots on the radius, both by 
plotting the mean values and standard deviations (error 
bars) of the systematic differences with respect to the 
input values of the simulations. The sum and the ratio 
of radii relative to the semimajor axis were used to check 
for systematic effects since these are the parameters that 
directly depend on the shape of the eclipses. A clear 
trend is seen in the case of Distribution 2, where the 
sum of the radii of the components (i.e., total size of 
the stars) is systematically overestimated by the fits by 
2-6%. For distributions less concentrated to the pole 
the deviations induced by spots seem to be random. In 
the case of the ratio of radii ratio, the differences found 
do not seem to be significant in any of the cases, and 
especially when considering the dispersion. The same 
lack of systematic trends is found for other fit parameters 
such as the inclination, effective temperature ratio, or 
light ratio. 

The parameters related to radial velocity curves, such 
as semimajor axis, mass ratio and systemic radial ve- 
locity, were also inspected for deviations by simulating 
the corresponding velocity curves. Fits using the WD 
code were performed by ignoring the presence of spots 
and yielded insignificant differences in the mass ratio and 
semi-major axis, in all cases below 1% and 0.5%, respec- 
tively. As a consequence, the absolute radius derived for 
each component of the eclipsing binary system is not bi- 
ased because of the effect of spots on the semi-major axis, 
which ultimately defines the scale of the system. 

The systematically larger radii found from the distribu- 
tions with polar spots have a relatively straightforward 



geometric interpretation. A projected star that has a 
polar cap loses its circular symmetry and its isophotes 
become elongated. This causes the eclipses to widen and 
the WD code, which assumes a Roche geometry (i.e., 
nearly spherical stars for these well-separated systems) 
finds a best fit with a larger stellar radius. In the case of 
spots more uniformly distributed in latitude, eclipses are 
not significantly affected, thus, systematic deviations ap- 
pear more random. We tested this by assuming systems 
with pure polar caps (i.e., circular and symmetric dark 
spots covering both caps) and carrying out fits without 
considering the presence of spots (since no modulations 
are seen). The results confirm the existence of the sys- 
tematic effect on the stellar radii found in the simulations 
with Distribution 2. Also, because of the much faster fits, 
we could test several input parameters (corresponding to 
the systems IM Vir, GU Boo, YY Gem, CU Cnc and 
CM Dra) that yielded comparable results, thus testing 
the systematic effect on radius also for different values of 
relevant binary parameters such as inclination, relative 
radius and temperature ratio. Using these simulations we 
also tested the influence of the photometric band by pro- 
ducing light curves in the BVRIJK filters. We found 
that the systematic effect on the radius is expectedly 
smaller in the near-IR bands because the spot contrast 
is lower. The radius differences range from ~ 3.5% in B 
to ^ 1.2% in K. The simultaneous fit of several bands 
still yields systematic radius differences, corresponding 
to some average of those resulting from fits to the indi- 
vidual filter light curves. 

The outcome of this battery of tests is that the pres- 
ence of starspots on the stellar surface could bias the 
determination of stellar radii. This is the case when 
spots are concentrated towards the poles. Doppler to- 
mography imaging has revealed a preeminence of polar 
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Fig. 6. — Fits to simulated liglit curves following the procedure explained in tlie text. Light curves are vertically shifted to be plotted 
with the same scale. Each block corresponds to a different spot distribution as indicated. Simulations with equal /3 value are arranged in 
columns, starting from /3 = 0.1 to 0.4 in steps of 0.1, except for Distribution 2, for which the last column corresponds to 13 = 0.3. 
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Fig. 7. — Differentials between input parameters to the simu- 
lations and those recovered from fits. The ratio and the sum of 
the radii are shown on the top and bottom panels, respectively, 
for different distributions: uniform (filled circles). Distribution 1 
(open squares) and Distribution 2 (filled triangles). These values 
are computed as the average of the results of fits on three simulated 
light curves for each distribution. Standard deviations are shown 
as error bars. 

spots on rapidly rotatin g low-mass stars, s uch as AG Dor, 
AB Dor, LO Peg (see iStrassmeieii [20091 for a review). 
YY Gem is the only EB system from our sample with 
published Doppler images ()Hatze£lll995f) , although no in- 
dication of a polar cap was found. However, polar spots 
are difficult to detect with Doppler tomography for near 
edge-on inclinations, such as the case of YY Gem. Spec- 
tropolarimetric analyses of low-m ass stars showing in- 
tense poloidal magnetic fields (see iDonati fc Landstreetl 
120091 for a review) may provide further evidence of po- 
lar spots. If polar caps are present on low-mass EBs, 
these could be responsible for up to ~ 6% (if (3 — 0.3) of 
the radius discrepancy between models and observations. 
Thus, stellar models should be compared with observed 
values corrected from the systematic effect due to polar 
spots. 

4. MODELLING ACTIVE STARS 

As already mentioned in Section [TJ IChabrier et al.l 
(|2007() introduced activity effects in low-mass stellar 
models by considering two parameters: the spot blocking 
factor (/?) and the modifica tion of the mixing len gth pa- 
rameter (a). The study of IChabrier et al] (|2007f) shows 
that the effects of the (3 and a parameters are degenerate, 
i.e., the properties of any given system can be reproduced 
by modifying any of the two or both, although for com- 
pletely convective stars, the effect of reducing the mixing 
length parameter is very small. 

CM Dra offers the possibility of discerning between the 
effects of activity on the convective efficiency (parame- 
terized by a) and on surface spots (parameterized by /?) . 
This is because the components of this system should 
be fully convective and their structure almost indepen- 
dent of the mixing length parameter a. Thus, we used 
this system to find the /? value that yields the best fit 
to the properties of the components of CM Dra. We 
compared the theoretical M — R relationships, interpo- 
lating between models with different /3. To do so, we 
assumed the presence of polar spots on the components 
(see Sect. 13. 2p . Different negative corrections to the radii 
of the components were used depending on the value of 



/3 and in consistency with our simulations in Fig. [71 
An iterative process was employed until reaching agree- 
ment between the resulting (3 value and the radius cor- 
rection from light curve systematics. After iterating, we 
find that a model with (3 = 0.17 ± 0.03 reproduces both 
components when their radii are downwards corrected by 
^ 3%. If spots are ~ 15% cooler than the photosphere, 
the resulting {3 translates into (36 ±6)% of the star sur- 
face covered by spots, a value that is in agreement with 
findings from Doppler imaging for other systems such as 
HK Aqr (Barnes & Collier Cameron 2001). Note, how- 
ever, that spot coverages from Doppler tomography rep- 
resent a lower limit because of the limited sensitivity of 
the technique to small-sized spots or low contrast tem- 
peratures. 

Abundant X-ray observations have revealed a clear re- 
lation between magnetic activity and rotation, indicating 
that rapidly rotating stars are more active than slower 
rotators up to a certai n limit where saturation is reached. 
iPizzolato et al.l ()2003D showed that M-dwarfs with rota- 
tion periods below ~ 10 days show saturation of their 
X-ray emission, therefore we could assume that binaries 
in our sample are all in this regime of saturated coronal 
activity. The physical mechanism of saturation is not 
well understood yet, but one of the feasible explanations 
could be that saturation is reached when the entire stellar 
corona is full of active regions. This may imply that the 
13 parameter is similar for stars that are in the saturated 
regime. 

At this point, we make two assumptions to proceed 
with the analysis: (1) that spots on active EBs are mostly 
polar and (2) that saturated systems have similar spot 
coverages (i.e., that /? is roughly independent of mass for 
very active systems). With these assumptions and fol- 
lowing the same procedure applied above for CM Dra, 
we find that the overall sample of EB systems is repro- 
duced within la of their error bars with (3 = 0.20 ± 0.04, 
i.e., (42 ± 8)% of spot coverage. 

The left panel of Figure [8| shows a comparison of mass 
and radius observations with models assuming different 
(3 values. The radii of the EB systems shown in the fig- 
ure were corrected considering two effects. First, a 3% 
systematic difference resulting from the light curve anal- 
ysis. Second, an age effect on the stellar radu to put all 
systems at a normalized age of 1 Gyr (ages in Table [l] ex- 
cept for NGC2204-S892 where no age is available and no 
correction was made) using the models of Baraffe et al.l 
(1998). The evolutionary radius offsets are always well 
below the 1% level. This allowed us to plot Figure [Sj and 
compare all EBs with a theoretical model with the same 
age. 

To serve as a reference, if we relax the polar spot as- 
sumption, i.e., we do not correct for the systematic effect 
of polar spots on light curves, the best model reproducing 
CM Dra gives f3 — 0.37 ± 0.04, which means more than 
75% of the surface of each component covered by spots. 
This seems to be a rather high value of spot coverage 
when compared with the results of Doppler imaging. If 
we do not consider the saturation hypothesis and assume 
that each system could be reproduced by different {3 fac- 
tors, we find in systems such as YY Gem or GU Boo over 
half of the surface must be spotted. Note that the spot 
coverage estimates depend on the spot contrast, which is 
set to ~ 500 K in this work. 
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Some remaining radius differences are still apparent 
in Figure [5] when considering a /3 value of 0.17. Systems 
such as YY Gem or GU Boo have observed radii well over 
Icr of the model predictions. It is interesting to note that 
these are the systems with the fastest rotational veloc- 
ity and perhaps be an indication of additional effects of 
rotation and/or magnetic activity not accounted for in 
the (3 analysis. The right panel of Figure [5] shows that 
for the more massive M-type EBs (< 0.6 M©) the addi- 
tional differences could be explained by reducing the mix- 
ing length parameter a. Figure |9] depicts the differences 
between the /3 = 0.17 model and observations (after cor- 
rection for polar spots and normalized to 1 Gyr) as a 
function of the rotational velocity of each star (computed 
assuming synchronization). Although not significant, a 
general trend of increasing differences for increasing ro- 
tational velocities seems to be present. The exception is 
NGC2204-S892 (two rightmost points on the plot), which 
is the shortest period binary. But note that no determi- 
nation of the age and metallicity of this EB system is 
available and this could be falsifying the analysis. The 
tentative trend may indicate that, while the spot effect 
(via the /3 parameter) is a clear contributor to explain- 



ing the differences between observation and theory, the 
loss of efficiency of convective energy transport (related 
with magnetic fields and thus possibly with rotation rate) 
may also be at play, explaining up to 4% of the radius 
difference. 

As stated in Section [21 absolute effective temperatures 
of EB components were not used to constrain model 
parameters because they are not determined indepen- 
dently of calibrations as masses and radii. Nevertheless 
we checked the consistency of the spot scenario with tem- 
peratures. The bottom panels of Figure [8] show the com- 
parison between models and observations in the M — Toff 
plane for different a and (3 values. The temperatures rep- 
resented on these plots have bee n also normahzed to an 
age of 1 Gyr using the models of iBaraffe et al.l (|1998f ) al- 
though differences amount for less than 10 K in all cases. 
The model with /3 — 0.17 also reproduces the M — Toff 
relationship of EBs within the errors, with the excep- 
tions of IM Vir and CU Cnc, whose temperatures may 
be affected by a circumstellar dust disk. These figures 
also show that both models with lower a and the model 
with /3 = 0.17 reproduces all systems but CM Dra, due 
to the negligible effect of mixing length on completely 
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Fig. 9. — Residual discrepancy between observations and the 
model with (3=0.17 that best fits the case of CM Dra, which is 
not affected by variations of the mixing length parameter a. 

convective stars. Thus, the conclusions from the M — R 
analysis are also consistent with the effective tempera- 
ture comparison. 

5. CONCLUSIONS 

After evaluating different scenarios to explain the dis- 
crepancies reported between observation and theory for 
low mass stars, a rather complex picture is now emerg- 
ing since no single effect can account for the full size of 
the differences. Our analysis shows that the ~ 5 — 10% 
radius discrepancy in M-type EBs is explained by com- 
bining three factors: (1) A systematic overestimation 
of the stellar radius in EB light curve analysis caused 
by the presence of polar spots with the subsequent loss 
of circular symmetry (amounting to about ~ 3%); (2) 
an increase of stellar radii to compensate for the loss of 
radiative efficiency because of starspots (reproduced by 
models considering a /3 parameter of 0.17 that explains 
2% of the radius difference); and (3) an increase of the 
radius caused by the lower convective efficiency in fast 
rotators with supposedly strong magnetic fields (repro- 
duced by models through changing the a parameter) and 
accounting for 0-4% of the radius difference. This sce- 
nario also reproduces the M — T^s distribution of EBs 
although effective temperatures are not as fundamentally 
determined as radii. 

Our overall analysis ensures self consistency by com- 
paring with the constraints set by the observations. Pre- 
dicted spot coverages and simulated light curves are 
found to reproduce the modulation amplitudes in the 
real light curves of the EB systems. Also the ^ 35% sur- 
face coverage is found to be compatible with results from 



Doppler tomography of active stars. Even relatively in- 
active systems such as some hosts of transiting exoplan- 
ets are found to have significant spo t coverages when 
applying transit mapping procedures (jPont et al.l 120071 : 
ISilva-Vahol[2008l) . 

It is important to note that our results are based on 
two main hypotheses: (1) That all of these active EB 
components have spots that preferentially occupy loca- 
tions close to the pole (with consequent effect on the 
light curve analysis) in line with results from Doppler 
imaging; and (2) that all EB components have similar 
surface spottedness caused by their saturated activity 
regime. The first hypothesis relaxes the classical dis- 
crepancy with models to a lower 2-7%, and the second 
hypothesis ensures that /3 will roughly be "universal" for 
these strongly active systems. In this framework, the 
system CM Dra is key in reaching the conclusions above. 
This is because it allows us to separate the effects of the 
two free parameters, namely /3 and a, since the compo- 
nents of this system are expected to be fully convective 
and little affected by changes in the mixing length pa- 
rameter. 

The analysis of low-mass stellar evolution models 
clearly shows that activity plays an important role on 
the structure of stars, and that when considering activity 
effects on stellar structure, models can reproduce the ob- 
served values of radii and effective temperatures of low- 
mass stars. Accurate masses and radii of new low-mass 
systems coming from missions such as CoRoT and Ke- 
pler should improve the statistics and better define the 
scenario proposed above, especially if comparable active 
and non-active EB systems are found. Doppler imag- 
ing of the EBs analyzed o r new methods of m apping 
stars using interferometry (|Monnier et al.l [20071 ) should 
also give indications on the spot distribution on the stel- 
lar surfaces. All these efforts are valuable to understand 
the vast majority of the stars that populate the Galaxy 
and, notably, have an impact in the area of exoplanet 
research, in which planet properties are directly linked 
to our knowledge of the host star. 
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